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1. Introduction 

A relativistic modified gravity (MOG) called Scalar- Tensor- Vector Gravity (STVG) 
desribes a self-consistent, stable gravity theory that contains Einstein's general rel- 
ativity in a well-defined limit P The theory has an extra degree of freedom, a vector 
field called a "phion" field whose curl is a skew symmetric field that couples to 
matter ("fifth force"). The spacetime geometry is described by a symmetric Ein- 
stein metric. An alternative relativistic gravity theory called Metric-Skew- Tensor 
Gravity (MSTG) has also been formulateoP in which the spacetime is described by 
a symmetric metric, and the extra degree of freedom is a skew symmetric second 
rank tensor field. Both of these theories yield the same weak field consequences for 
physical systems. 

The classical STVG theory allows the gravitational coupling "constant" G and 
the coupling of the phion field and its effective mass to vary with space and time 
as scalar fields. 

A MOG should explain the following physical phenomena: 
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(1) Galaxy rotation curve data; 

(2) Mass profiles of x-ray clusters; 

(3) Gravitational lensing data for galaxies and clusters of galaxies; 

(4) The cosmic microwave background (CMB) including the acoustical oscillation 
power spectrum data; 

(5) The formation of proto-galaxies in the early universe and the growth of galaxies; 

(6) N-body simulations of galaxy surveys; 

(7) The accelerating expansion of the universe. 

We seek a unified description of solar system, astrophysical and large-scale cos- 
mological data without exotic non-baryonic dark matter. Dark matter in the form 
of particles has until now not been discovered in spite of large-scale experimental 
efforts.^ The accelerating expansion of the universe should be explained by the 
MOG theory without postulating a cosmological constant. 



2. Action and Field Equations 

Our MOG action takes the fornffl 



where 



and 
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(5) 
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We have chosen units with c = 1, V M denotes the covariant derivative with 
respect to the metric g^ u . We adopt the metric signature i]^ = diag(l, — 1, — 1, — 1) 
where rj^ is the Minkowski spacetime metric, and R — g^ 1 'Ruv where R^ v is the 
symmetric Ricci tensor. Moreover, V(<j)) denotes a potential for the vector field </> M , 
while V(G),V(ui) and V(fi) denote the three potentials associated with the three 
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scalar fields G(x),u)(x) and fi(x), respectively. The field u>(x) is dimensionless and 
A denotes the cosmological constant. Moreover, 

B^ v = d^ v ~ d u ^,j,. (8) 

The field equations and the test particle equations of motion are derived in Ref. [TJ 
The action for the field is of the Maxwell-Proca form for a massive vector 
field cf)^. It can be proved that this MOG possesses a stable vacuum and the Hamil- 
tonian is bounded from below. Even though the action is not gauge invariant, it can 
be shown that the longitudinal mode 4>$ (where </) M = (</>o, <f>i) (i = 1, 2, 3)) does not 
propagate and the theory is free of ghosts. Similar arguments apply to the MSTG 
theoryll] 



3. Modified Newtonian Acceleration Law and Galaxy Dynamics 

The modified acceleration law can be written aiffl 

G{r)M 



where 



G(r) = G 



N 



a(r) 



I Mo 



M 



— l-exp(-r/r ) 1 + — 



''o 



(9) 



(10) 



is an effective expression for the variation of G with respect to r, and Gjv denotes 
Newton's gravitational constant. A good fit to a large number of galaxies has been 
achieved with the parameters^ 



M = 9.60 x 10 11 M Q , r = 13.92 kpc = 4.30 x 10^ cm 



22 



(11) 



In the fitting of the galaxy rotation curves for both LSB and HSB galaxies, using 
photometric data to determine the mass distribution M. (r) pi only the mass-to-light 
ratio (M/L) is employed, once the values of M and are fixed universally for all 
LSB and HSB galaxies. Dwarf galaxies are also fitted with the parameters^ 



M = 2.40 x 10 U M Q , r = 6.96 kpc = 2.15 x 10 22 cm. 



(12) 



By choosing universal values for the parameters Goo = Gjv(1 + y/Mo/M), (Mo) dust 
and (rp)ciust, we are able to obtain satisfactory fits to a large sample of X-ray cluster 
dataP 



4. Solar System and Binary Pulsar 

Let us assume that we are in a distance scale regime for which the fields G, u> and 
/i take their approximate renormalized constant values: 

G~G (1 + Z), w~uoA, m~W>5, (13) 

a For a detailed discussion of possible instabilities and pathological behavior of vector-gravity the- 
ories, see Ref. [4] 
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where Go,wo & n d fJ-o denote the "bare" values of G,lo and /i, respectively, and Z, A 
and B are the associated renormalization constants. We obtain from the equations 
of motion of a test particle the orbital equation (we reinsert the speed of light c 

-^^u 2 . (14) 



d 2 u GM K 

d4 : 



2 1 u = 7Tn ~^12 c M-r/r ) 



1+1- 



where u = 1/r, K = Gn^/MMq and J denotes the orbital angular momentum. 
Using the large r weak field approximation, we obtain the orbit equation for r<ro: 

d 2u »r „ GM 2 

= N + 3^-u 2 , (15) 



:,2 



where Jm denotes the Newtonian value of J and 

, T GM K 

N =^J2-^JT- ( 16 ) 

We can solve Eq. (|15p by perturbation theory and find for the perihelion advance 
of a planetary orbit 

Ac = ^L(GM G -K ), (17) 

C Li 

where Jn — (GMqL/c 2 ) 1 / 2 , L = a(l — e 2 ) and a and e denote the semimajor axis 
and the eccentricity of the planetary orbit, respectively. 

For the solar system r -C tq and from the running of the effective gravitational 
coupling constant, G = G(r), we have G ~ Gat within the experimental errors for 
the measurement of Newton's constant Gn- We choose for the solar system 

1.5km (18) 

cr 

and use G = G^r to obtain from (|17p a perihelion advance of Mercury in agreement 
with GR. The bound (fT8|) requires that the coupling constant oj varies with distance 
in such a way that it is sufficiently small in the solar system regime and determines 
a value for Mq that is in accord with the bound (fT8|) . 

For terrestrial experiments and orbits of satellites, we have also that G ^ Gn and 
for Kq sufficiently small, we then achieve agreement with all gravitational terrestrial 
experiments including Eotvos free-fall experiments and "fifth force" experiments. 

For the binary pulsar PSR 1913+16 the formula (fT7)l can be adapted to the 
periastron shift of a binary system. Combining this with the STVG gravitational 
wave radiation formula, which will approximate closely the GR formula, we can 
obtain agreement with the observations for the binary pulsar. The mean orbital 
radius for the binary pulsar is equal to the projected semi-major axis of the binary, 
(r) n = 7 x 10 10 cm, and we choose (r) n "C tq. Thus, for G = Gat within the exper- 
imental errors, we obtain agreement with the binary pulsar data for the periastron 
shift when 

^«4.2km. (19) 
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For a massless photon we have 

d 2 u 



GM 2 
u = 3 — 7T-u 



(20) 



For the solar system using G ~ Gjv within the experimental errors gives the light 
deflection: 

4G N M Q 



in agreement with GR. 



A = 



c 2 R, 



(21) 



5. Pioneer Anomaly 

The radio tracking data from the Pioneer 10/11 spacecraft during their travel 
to the outer parts of the solar system have revealed an anomalous acceleration. 
The Doppler data obtained at distances r from the Sun between 20 and 70 astro- 
nomical units (AU) showed the anomaly as a deviation from Newton's and Ein- 
stein's gravitational theories. The anomaly is observed in the Doppler residuals 
data, as the differences of the observed Doppler velocity from the modelled Doppler 
velocity, and can be represented as an anomalous acceleration directed towards 
the Sun, with a n appro ximately constant amplitude over the range of distance, 
20 au < r < 70 a- 



tnapprox 
u 7l8l9ll01 



a P = (8.74 ± 1.33) x 10~ 8 cms~ 



(22) 



After a determined attempt to account for all known sources of systematic errors, 
the conclusion has been reached that the anomalous acceleration towards t he Sun 
could be a real physical effect that requires a physical explanationP^M It f> 
We can rewrite the acceleration in the form 

G N M 



a{r) 



1 + a(r) 



1 - exp(-r/A(r)) 1 + 



A(r) 



(23) 



We postulate a gravitational solution that the Pioneer 10/11 anomaly is caused 
by the difference between the running of G(r) and the Newtonian value, Gn- So 
the Pioneer anomalous acceleration directed towards the center of the Sun is given 

by 



ap 



SG(r)M Q 



where 



SG(r) = G N a(r) 



1 - exp(-r/A(r)) 1 + 



A(r) 



(24) 



(25) 



b It is possible that a heat transfer mechanism from the spacecraft transponders could produce a 
non-gravitational explanation for the anomaly. 
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Lacking at present a solution for the variations of a(r) and A(r) in the solar system, 
we adopt the following parametric representations of the "running" of a(r) and 
A(r): 

a{r) = a<30 (l - exp(-r/r)) b/2 , (26) 

A(r) = A^Cl-exp^r/f))- 6 . (27) 

Here, f is a non-running distance scale parameter and b is a constant. 

In Ref. [Til a best fit to the acceleration data extracted from Figure 4 of Ref. 10 
was obtained using a nonlinear least-squares fitting routine including estimated 
errors from the Doppler shift observations". The best fit parameters are 

aoo = (1.00 ±0.02) x 10~ 3 , 
Aoo = 47 ± 1 au, 
f = 4.6 ± 0.2 au, 

b = 4.0. (28) 

The small uncertainties in the best fit parameters are due to the remarkably low 
variance of residuals corresponding to a reduced x 2 per degree of freedom of 0.42 
signalling a good fit. An important result obtained from our fit to the anomalous 
acceleration data is that the anomalous acceleration kicks-in at the orbit of Saturn. 

Fifth force experimental bounds plotted for log 10 a versus log 10 A are shown in 
Fig. 1 of Ref. fT2l for fixed values of a and A. The updated 2003 observational data 
for the bounds o bta ined from the planetary ephemerides is extrapolated to r = 
10 15 m = 6, 685 arPl However, this extrapolation is based on using fixed universal 
values for the parameters a and A. Since known reliable data from the ephemerides 
of the outer planets ends with the data for Pluto at a distance from the Sun, 
r = 39.52 au = 5.91 x 10 12 to, we could claim that for our range of values 47 au < 
A(r) < oo, we predict a(r) and A(r) values consistent with the un- extrapolated fifth 
force bounds. 

A consequence of a variation of G and GMq for the solar system is a modification 
of Kepler's third law: 

a PL =G(a PL )M G (^^j , (29) 

where T P l is the planetary sidereal orbital period and ap^ is the physically mea- 
sured semi-major axis of the planetary orbit. For given values of apL and Tpl, (|29j) 
can be used to determine G(r)M Q . 

For several planets such as Mercury, Venus, Mars and Jupiter there are plan- 
etary ranging data, spacecraft tracking data and radiotechnical flyby observations 
available , and i t is possible to measure ap^ directly. For a distance varying GAIq 
wo deriving 

"G(a PL )M l 1/3 



(^1=1 

\apL, 



Vpl 



G(a e )M 



(30) 
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Here, it is assumed that GMq varies with distance such that rjpL can be treated as 
a constant for the orbit of a planet. We obtain 

[ G(a P£ ) l 1/3 

VPL = -jt, r- - 1. (31) 

L G{a @ ) . 

The results for At/pl due to the uncertainty in the planetary ephemerides 
are presented in Ref. |TT] for the nine planets and are consistent with the solar 
ephemerides. 

The validity of the bounds on a possible fifth force obtained from the ephemerides 
of the outer planets Uranus, Neptune and Pluto are critical in the exclusion of a 
parameter space for our fits to the Pioneer anomaly acceleration. Beyond the outer 
planets, the theoretical prediction for rj(r) approaches an asymptotic value: 

^ = lira r)(r) = 3.34 x 10" 4 . (32) 

r — >oc 

We see that the variations ("running") of a(r) and A(r) with distance play an 
important role in interpreting the data for the fifth force bounds. This is in contrast 
to the standard non-modified Yukawa correction to the Newtonian force law with 
fixed universal values of a and A and for the range of values < A < oo, for which the 
equivalence principle and lunar laser ranging and radar ranging data to planetary 
probes exclude t he possib ility of a gravitational and fifth force explanation for the 
Pioneer anomaly! I I I 

A study of the Shapiro time delay prediction in our MOG is found to be con- 
sistent with time delay observations and pre dicts a measurable deviation from GR 
for the outer planets Neptune and PlutoM^ 

6. Gravitational Lensing 

The bending angle of a light ray as it passes near a massive system along an ap- 
proximately straight path is given to lowest order in v 2 /c 2 by 

2 ' \a^\dz, (33) 



c 2 



where _L denotes the perpendicular component to the ray's direction, and dz is the 
element of length along the ray and a denotes the acceleration. 
From (|20|) . we obtain the light deflection 

. 4GM AG N M 

A = ^ = ^iT' (34) 

where 

M = M(l + S J^y (35) 
The value of M follows from (fTU)) for clusters as r > r and 

G(r)^G oc = G N (l + J^-). (36) 
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We choose for a cluster Mq = 3.6 x 10 15 M© and a cluster mass M c i us t ~ 10 14 Mq, 
and obtain 

/ / W„ \ 

~ 6. (37) 

clust 

We see that M ~ 7M and we can explain the increase in the light bending without 
exotic dark matter. 
For r > ro we get 

G N M 

a{r) = — . (38) 

We expect to obtain from this result a satisfactory description of lensing phenomena 
using Eq.i 



7. Modified Friedmann Equations in Cosmology 

We shall base our results for the cosmic microwave background (CMB) power spec- 
trum on our MOG without a second compo nent of cold dark matter (CDM). Our 
description of the accelerating universe ^^^ is based on Aq in Eq. ([62)) derived from 
our varying gravitational constant!^ 

We adopt a homogeneous and isotropic Friedmann-Lemaitre- Robertson- Walker 
(FLRW) background geometry with the line element 

ds 2 = dt 2 -aW^^+rWY (39) 

where d£l 2 = d6 2 + sin 2 Odtfi 2 and k = 0, — 1, +1 for a spatially flat, open and closed 
universe, respectively. Due to the symmetry of the FLRW background spacetime, 
we have 4>q = <p ^ 0, <pi = and = 0. 

We define the energy-momentum tensor for a perfect fluid by 

= {p+p)vPu v -psT, (40) 

where u M = dx^ /ds is the 4-velocity of a fluid clement and g tlv u tl u 1 ' — 1. Moreover, 
we have 

P = Pm+ P4> + PS, P=Pm+P<t>+PS, (41) 

where pi and pi denote the components of density and pressure associated with the 
matter, the <\>^ field and the scalar fields G, lo and ^.respectively. 
The modified Friedmann equations take the forrrW 

«»(t) k _ 8nG(t)p(t) A 
a*(t) + a2(t)~ 3 +'W+ 3 ' 

^ = -^W*)+3p(*)] + M*) + ^ (43) 
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where a = da/dt and 
From (l42l we obtain 



Pa3 = 8^G a + k ~ ° 2/ " ■ (46) 

This leads by differentiation with respect to t to the expression: 

d In a , , _ . 
p + 3— (p+p)+X = 0, (47) 

where 

3a 2 

l=—(2af + af-2dh). (48) 

An approximate solution to the field equations for the variation of G in Ref. Q] 
in the background FLRW spacetime is given by 

g + 3Hg + V\g) = ±G N g 2 (p-3p +I ^r-^j, (49) 

where Q(t) = G{€)/Gn and H — a/a. A solution for Q in terms of a given potential 
V{g) and for given values of p, p and A can be obtained from ([49J) 122J 

The solution for g must satisfy a constraint at the time of big bang 
nucleosynthesis!^ The number of relativistic degrees of freedom is very sensitive to 
the cosmic expansion rate at 1 MeV. This can be used to constrain the time depen- 
dence of G. Measurements of the A He mass fraction and the deuterium abundance 
at 1 MeV lead to the constraint G{t) ~ Gjv- We impose the condition Q(t) — * 1 as 
t — ► tg bn where ts bn denotes the time of the big bang nucleosynthesis. Moreover, 
locally in the solar system we must satisfy the observational bound from the Cassini 
spacecraft measurements^! 

\G/G\ < lO-^yr -1 . (50) 
We shall now impose the approximate conditions at the epoch of recombination: 



2af + af ~2d/i, (51) 

m < 41 



We find from (|45| and (|52jl that / ~ h, and from the condition ijSlj) we obtain 

/-^~0, (53) 
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where 

A C ^§. ,54, 

By setting the cosmological constant A = 0, we get the generalized Friedmann 
equations 

a 2 k 8ttGp 



a 2 a 2 3 
a 47rG 



Ac, (55) 



„ (p + 3 P ) + A G . (56) 
a 3 

We now have from (|4"T)) , l|48p and (fSTj) at the epoch of recombination X ~ and 

p + 3^(p+p)~0. (57) 

We adopt the equation of state: p(t) = wp(t) and derive from (|5T[) the approximate 
solution for p(i): 

p(*W(*0(^) , (58) 

where a/a — 1/(1 + z) and z denotes the red shift. For the matter and radiation 
densities p m and p r , we have w = and u> = 1/3, respectively. This gives 

p m (t) ~ Pm {t )(l + z) 3 , p r (t) ~ p r (t )(l + z) 4 . (59) 

Let us expand G(£) in a power series 

= G cS (U) + (t- t r )G{t r ) + {t- t r ) 2 G{t r ) + .... (60) 

where t ~ i r is the time of recombination and G e ff(£ r ) = Gjv(1 + Z) = const. We 
write the generalized Friedmann equation for flat space, k — 0, in the approximate 
form 

H 2 = 8nGcSPm +A G , (61) 
3 

where 

A G = > (62) 



and Aq ~ 0. It follows from (|6T|) that for a spatially flat universe: 

n m + n G = i, (63) 

where 

87rG ffPm A G 

We shall postulate that the matter density p m is dominated by the baryon density, 
Pm ~ Pb, and we have 

n TO ~ Obeff) (65) 



4, 2008 13:9 WSPC/INSTRUCTION FILE nasavirgproc5 



A modified gravity and its consequences for the solar system, astrophysics and cosmology 11 

where 

d>irG cS p b 

S = 3H 2 ' ^ ' 

Thus, we assume that the baryon-photon fluid dominates matter before recombina- 
tion and at the surface of last scattering without a cold dark matter fluid component. 

From the current value: Ho = 7.5 x 10~ n yr" 1 and (|62l) and (|64|) . we obtain for 
n G ~ 0.7: 

\G/G\ -5 x 10 _11 yr _1 , (67) 

valid at cosmological scales for red shifts z > 0.1. In the local solar system and for 
the binary pulsar PSR 1913+16 for z ~ 0, the experimental bound is 

\G/G\ < 5 x l(T 12 yr _1 - (68) 

We can explain the accelerated expansion of the universe deduced from supernovae 
measurements in the range 0.1 < z < 1.7 using the cosmologically scaled value of 
G/G in l|67p with Einstein's cosmological constant A = 0. 



8. Acoustical Peaks in the CMB Power Spectrum 

Mukhanov^S has obtained an analytical solution to the amplitude of fluctuations 
in the CMB power spectrum for I 3> 1: 



l{l + l)C t ~ -(O + N) 



(69) 



Here, O denotes the oscillating part of the spectrum, while the non-oscillating con- 
tribution can be written as the sum of three parts 



^ = ^1+^2+ N 3 . 

The oscillating contributions can be calculated from the formula 



0~J — 

Thl 



Ax cos [ lr p + — ) + A-2 cos 2lr p + — 



exp(-(;/Z s ) 2 ), 



(70) 



(71) 



where r^ and r p are parameters that determine predominantly the heights and 
positions of the peaks, respectively. The A\ and A-i are constant coefficients given 
in the range 100 < I < 1200 for Q m ~ QbeB by 

X(P -0.78) 2 -4.3) 



A 1 ~0.1£- 



(1 + 1/4 



lf)l 2 



A 2 ~ 0.14 



(0.5 + 0.36P) 2 
(1+0 1/2 



where 



V = In 



II 

200(f! hcff )i/2;' 



(72) 
(73) 

(74) 
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and / is given by the ratio 



i 

J/2 



9.C 



bell 



1/6 



dx 



o (sinhx) 2 /^ Z V2 



(75) 



Here, rj x and z x denote a conformal time r\ — r\ x and a redshift in the range rj > 
Vx > Vr when radiation can be neglected and y = sinh^ 1 (Sla/Qbes) 1 / 2 ■ To determine 
Vx/vo, we use the exact solution for a flat dust-dominated universe with a constant 
A G : 

a(t) =a (smh(^\H t\ , (76) 

where oq and Hq denote the present values of a and the Hubble parameter H. 

The If and l s in (f?2"j) denote the finite thickness and Silk damping scales, re- 
spectively, given by 



,2 _ 1 fm 

f ~ 2a 2 {r, r 



H = 



1 



no 



where 

a ~ 1.49 x 10~ 2 



1 



1 



■°cq 



2(a 2 + l/(k DVr ) 2 ) \r) r 

-1/2- 



2 

5 Jo 



-1/2 



(77) 



(78) 



and r 7 is the photon mean-free time 
A numerical fitting formula givei 



In 



I 



2oo(n° e | 9 ) 



— / dr]c s (r/). 



Moreover, 



3 ( Pb 



^ 3c 2 1 4Vp 7 



where c s {rf) is the speed of sound: 



c s {v) 



1 

V3 



a(r) r ) 



1-1/2 



(79) 



(80) 



(81) 



We note that £ does not depend on the value of G e ff- For the matter-radiation 
universe: 

2 



a(rj) = a 



where for radiation-matter equality z = z eq : 

2 



— ~ I — ) +2[ ^ 



(82) 



(83) 



and rjeq = ^*(\/2 — 1) follows from a = a{rj eq ). 
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For the non-oscillating parts, we have 

^ (1 ^/-S^ J)2 ^-"/-.)'). 

0.033 7>-0.5(7/U a55 + 2.2) 2 . ,, ,, n2n 

^ ~ (TTO^ i + 2 (W ^-^-) a )- (86) 

Mukhanov's formula^] for the oscillating spectrum is given by 

where we have normalized the power spectrum by using for a flat spectrum with a 
constant amplitude B: 

9B 



[i(i+l)Q]towJ = 



IOOtt 

We adopt the parameters 

n bN ~om, f> bcff ~o.3, f> G ~o.7, £~o.6, ( 8 g) 

and 

r, l= 0.03, r p = 0.01 l } ~ 1580, Z s - 1100, (90) 

where fi^jv = SnGNPh/SH 2 . 

The fluctuation spectrum determined by Mukhanov's analytical formula is dis- 
played in Fig. 1 for the choice of cosmological parameters given in ([89)) and (|90| . 



The role played by CDM in the standard scenario is replaced in the modified 
gravity theory by the significant deepening of the gravitational potential well by 
the effective gravitational constant, G e s ~ 7Gn, that traps the non-relativistic 
baryons before recombination. The deepening of the gravitational well reduces the 
baryon dissipation due to the photon coupling pressure and the third and higher 
peaks in the acoustic oscillation spectrum are not erased by finite thickness and 
baryon drag effects. The effective baryon density flbeff = (1 + Z)Q&jy ~ 7&bN ~ 0.3 
dominates the fluid before recombination, and we fit the acoustical power spectrum 
data without a cold dark matter fluid component. For t < tdoc, where td oc denotes 
the time of matter-radiation decoupling, luminous baryons and photons are tightly 
coupled and for photons the dominant collision mechanism is scattering by non- 
relativistic electrons due to Thompson scattering. It follows that luminous baryons 
are dragged along with photons and perturbations at wavelength X w < £ s will be 
partly erased where l s is the proper Silk length given by £ s ~ 3.5 Mpcf2 6e ]j !^ We 
have i s ~ 6Mpc for Q bcS ~ 0.3 compared to i s ~ 18Mpc for Q bN ~ 0.04. The Silk 
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Fig. 1. The solid line shows the result of the calculation of the power spectrum acoustical os- 
cillations: C(l), and the os correspond to the WMAP, Archeops and Boomerang data in units 
liK 2 X 10~ 3 as presented in Refs. 27-30. 

mass is reduced by more than an order of magnitud^. Thus, sufficient baryonic 
perturbations should survive before t ~ td cc to explain the power spectrum without 
collisionless dark matter. 

Our predictions for the CMB power spectrum for large angular scales corre- 
sponding to I < 100 will involve the integrated Sachs- Wolfe contributions obtained 
from the modified gravitational potential. 



9. Conclusions 

m 

We have demonstrated that a modified gravity theory J can lead to a satisfactory 
fit to the galaxy rotation curve data, mass profiles of x-ray cluster data, the so- 
lar system and the binary pulsar PSR 1913+16 data. Moreover, we can provide 
an explanation for the Pioneer 10/11 anomalous acceleration data, given that the 
anomaly is caused by gravity. We can fit satisfactorily the acoustical oscillation 
spectrum obtained in the cosmic microwave background data by employing the 
analytical formula for the fluctuation spectrum derived by MukhanovM^I 

Aq obtained from the varying gravitational constant in our MOG replaces the 
standard cosmological constant A in the concordance model. Thus, the accelerating 
expansion of the universe is obtained from the MOG scenario. 

An important problem to investigate is whether an N-body simulation calcula- 
tion based on our MOG scenario can predict the observed large scale galaxy surveys. 
The formation of proto-galaxy structure before and after the epoch of recombination 
and the growth of galaxies and clusters of galaxies at later times in the expansion 
of the universe has to be explained. 



c Notc that there will be a fraction of dark baryonic matter before decoupling. 
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We have succeeded in fitting in a unified picture a large amount of data over 16 
orders of magnitude in distance scale from Earth to the surface of last scattering 
some 13.7 billion years ago, using our modified gravitational theory without exotic 
dark matter. The data fitting ranges over four distance scales: the solar system, 
galaxies, clusters of galaxies and the CMB power spectrum data at the surface of 
last scattering. 
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